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8p22-23-rs2254546 was firstly discovered to be associated with Kawasaki disease (KD) susceptibility by a 
genome-wide association study. However, only one Chinese replication study has been performed so far. To 
verify this association in another Chinese population, a hospital-based case-control study in Zhejiang 
province was conducted followed by an integrated meta-analysis, comprising five case-control studies of 
1958 cases, 5615 controls and four transmission disequilibrium tests of 503 trios. In our case-control study, 
significant associations were observed between GG genotype or GG/GA genotypes of rs2254546 and 
increased KD risk (OR - 1.86, 95% CI - 1.01-3.41, P - 0.045; OR - 1.83, 95% CI - 1.01-3.33, P - 0.048), 
compared with AA genotype; however, no significant association was found in allelic model (OR — 1.20, 
95% CI — 0.96-1.50, P — 0.117). The meta-analysis further revealed that the G allele was significantly 
associated with the increased KD risk without evidence of heterogeneity (OR — 1.55, 95% CI — 1.42-1.70, P 
< 0.001). In conclusion, rs2254546 polymorphism might significantly contribute to the risk of KD. 

Kawasaki disease (KD; OMIM, 611775), also called Mucocutaneous Lymph Node Syndrome, is an acute, 
self-limited vasculitis that predominately affects infants and young children' Although it was first 
described in Japan in 1967 by the Japanese pediatrician Tomisaku Kawasaki'*, the etiology of KD is stiU 
not completely clear. The major clinical manifestations of KD include prolonged fever, bilateral non-purulent 
conjunctivitis, diffuse mucosal inflammation, polymorphous skin rashes, peripheral extremity changes, and 
cervical lymphadenopathy^. About 15-25% of untreated cases suffer coronary artery lesions, which is the most 
serious complication of KD, and makes KD the leading cause of acquired heart disease in children of developed 
countries'". 

Genetic determinants have been suggested to contribute to KD susceptibility. Firstly, the risk of KD in siblings 
of affected children is ten times higher as compared with the general population'. Secondly, it is twice as high in 
children born to parents with a history of KD compared with the general population". Thirdly, populations in 
Asian countries have higher incidence rates of KD than those in Western countries: Japan has the highest annual 
incidence rate'-", followed by Korea'" and Taiwan", regardless of their residences'^. Accumulating clinical and 
epidemiological findings suggest that infectious agents trigger inflammatory responses, stirring up host immune 
dysregulation in part of genetically predisposed individuals" '''. 

Genome-wide association studies (GWASs) have revolutionized the identification of genomic regions assoc- 
iated with complex diseases'^, and so do with KD. Until now, a total of sfac published GWASs have identified 
multiple novel single nucleotide polymorphisms (SNPs) associated with KD or coronary artery lesions in 
Japanese, Europeans, Koreans, Taiwanese and so on"'"^'. Among these SNPs, rs2254546, located at an intergenic 
region between family with sequence similarity 167, member A {FAM167A, also known as C8orfl3) and B 
lymphoid tyrosine kinase (BLK) at 8p22-23, was firstly discovered to be associated with KD by Onouchi 
et al.^°, who firstly performed a GWAS in a large sample size and then validated this locus in two independent 
case-control panels of Japanese and transmission disequilibrium tests (TDTs) of other ancestries. Notably, the 
TDT based on family is considerably advantageous as less confounding caused by population admixture and is 
more important than case-control study in genetic association analysis^^. 
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The results of KD research on Chinese are very important because 
of the high prevalence, but as far as we know, there is only one study 
which has replicated the association between rs2254546 and KD risk 
on a population from Sichuan province of Southwest China^"*. 
Therefore, we proposed to verify this association on another 
Chinese population from Zhejiang province of Southeast China. 
Subsequently, a meta-analysis, combining the current and previously 
published studies, and applying the method described by Kazeem 
and FarraU^**, was performed to integrate the results from case-con- 
trol studies and TDTs, and to provide more convincing evidence for 
this association. 

Results 

Results of our current case-control study. Population characteris- 
tics. A total of 428 cases and 478 controls participated in our current 
study. There were 263 (61.45%) males in cases and 294 (61.51%) 
males in controls, and no statistically significant difference was 
observed between these two groups in the distribution of gender 
(Pearson = 0.000, P = 0.986). 

Association analysis. Among all the participants, 892 samples were 
genotyped successfully, so that the genotyping call rate of rs2254546 
was 98.45%. The distribution of genotypes in controls was in agree- 
ment with Hardy- Weinberg equilibrium (HWE, P = 0.062) while no 
significant difference was observed in genotype frequencies between 
cases and controls (x^ = 4.08, P = 0.130). As shown in Table 1, under 
the multivariate logistic regression model adjusted for gender, indi- 
viduals with GG genotype showed a significant association with the 
increased risk of KD, compared with those carrying AA genotype 
(odds ratio (OR) = 1.86, 95% confidence interval (95%CI) = 1.01- 
3.41, P = 0.045). The same result was found under the dominant 
model (OR = 1.83, 95% CI = 1.01-3.33, P = 0.048). However, we 
failed to find the significant association in allelic model (OR = 1.20, 
95% CI = 0.96-1.50, P = 0.117). 

Results of meta-analysis. Study characteristics. After comprehensive 
searching, two published articles met the inclusion criteria, one 
GWAS and one replication study respectively^"'^^. Thus, a total of 
eight previous reports plus our current study, comprising five case- 
control studies of 1958 cases, 5615 controls and four TDTs of 503 
trios were finally included in this meta-analysis. Of the nine studies, 
six studies were conducted in Asians, one in Europeans, one in 
Hispanics and the other one in multi-ethnic population. Table 2 
shows the characteristics of the included studies. 

Combined results of case-control studies and TDTs. No significant 
evidence of heterogeneity was observed among the five case-control 
studies and four TDTs (Q = 13.00, -Phcterogeneity = 0.112); thus, the 
fixed-effects model was employed to assess the pooled association 
between rs2254546 and KD susceptibility. In the overall meta-ana- 
lysis, the G allele conferred an increased risk of KD compared with 
the A allele (OR = 1.55, 95% CI = 1.42-1.70, P < 0.001; Figure 1). 



Stratified analysis. In consideration of our strategy of meta-analysis 
by integrating case-control studies and TDTs, the data were firstly 
stratified by study design. As shown in Table 3, significant 
heterogeneity was detected in case-control subgroup but not in 

TDT subgroup (Q = 10.84, i^heterogenelty ~ 

0.028; Q = 0.83, 

J'heterogencity = 0.843). MeanwhUe, significant associations were 
found between the polymorphism and increased KD risk in both 
case-control and TDT subgroups (OR = 1.57, 95% CI = 1.34- 
1.85, P < 0.001; OR = 1.35, 95% CI = 1.05-1.74, P = 0.018), 
under random-effects model and fixed-effects model respectively. 

We further stratified the data by ethnicity. Because there was less 
than three studies conducted in Europeans or Hispanics, we per- 
formed a meta-analysis in Asian population only. Similarly, we 
found a significant association between the G allele of rs2254546 
and increased KD risk under random-effects model due to the het- 
erogeneity (OR = 1.58, 95% CI = 1.44-1.74, P < 0.001, Q = 10.87, 

f'heterogeneity = 0.054; Table 4). 

Sensitivity analysis. To evaluate the influence of each study on the 
overall estimate, sensitivity analysis was performed. A series of 
pooled ORs with 95% CIs produced similarly before and after 
eliminating each study at a time, suggesting that our results were 
robust (Table 5). 

Cumulative meta-analysis. Cumulative meta-analysis was con- 
ducted via adding one study at a time in chronological order. 
Although the cumulative effects tended to be null significant in 
time sequence, they still kept exceeding the threshold 1.0. 
Meanwhile, the 95% CIs became increasingly narrower with the 
accumulation of more data, indicating a progressive improvement 
in the precision of the estimates (Figure 2). 

Publication bias. Publication bias was firstly assessed by visual 
inspection of the funnel plot. The dots were nearly symmetrically 
distributed, suggesting little evidence of publication bias (Figure 3). 
The result of regression test for funnel plot asymmetry further 
indicated no pubhcation bias in this meta-analysis (P = 0.815). 

Discussion 

In our current case-control study, individuals carrying GG genotype 
or GG/GA genotypes of rs2254546 showed significant associations 
with KD susceptibility, compared with AA genotype. Moreover, the 
following meta-analysis, integrating five case-control studies and 
four TDTs, provided better evidence of the significant relationship 
between the polymorphism and increased KD risk, which was fur- 
ther confirmed by cumulative analysis in chronological order. Since 
the discrepancy in genetic backgrounds has been implied to be one 
possible reason for the distinction of diseases or traits in different 
ethnicities^\ and the prevalence of KD is overwhelmingly higher in 
Asians, it was deemed meaningful to assess the association of 
rs2254546 with KD risk in Asian ethnicity alone. Besides, sensitivity 
analysis indicated the stability of the results and little publication bias 
existed. 



Table 1 | Association between rs2254546 and KD 
Genotype Cases N (%) 


risk in the current study 
Controls N {%) 


OR 


95%CI 


Pi 


AA 


1 7 (4.07) 


34 (7.17) 


Ref. 


Ref. 




AG 


138 (33.01) 


156 (32.91) 


1.77 


0.95-3.31 


0.073 


GG 


263 (62.92) 


284 (59.92) 


1.86 


1.01-3.41 


0.045 


G/A allele 






1.20 


0.96-1.50 


0.117 


Dominant model 






1.83 


1.01-3.33 


0.048 


Recessive model 






1.14 


0.87-1.49 


0.347 


Additive model 






1.19 


0.95-1.48 


0.123 


^The Pvaiues were adjusted for gender. 
The positive results were in bold. 
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Table 2 | The characteristics of included studies in meta-analysis 



First author 


Year 


Country 


Ethnicity 


Study design 


Study type 


Genotyping method 


Sample size 


Onouchi 


201 2 


Japan 


Asian 


CC 


GWAS 


lllumina Human Hap550v3 BeadChip 


428/3379 




201 2 






CC 


IxCkJI II^U 1 l(J1 1 




470/378 


Onouchi 


2012 


Japan 


Asian 


CC 


Replication 


Invader Assay 


284/569 


Onouchi 


2012 


USA 


European 


TDT 


Replication 


Invader Assay 


248 


Onouchi 


2012 


USA 


Asian 


TDT 


Replication 


Invader Assay 


72 


Onouchi 


2012 


USA 


Hispanic 


TDT 


Replication 


Invader Assay 


94 


Onouchi 


2012 


USA 


Mixed 


TDT 


Replication 


Invader Assay 


89 


Yan 


2013 


China 


Asian 


CC 


Replication 


TaqMan SNP Genotyping Assays 


358/815 


Current study 


2013 


China 


Asian 


CC 


Replication 


TaqMan SNP Genotyping Assays 


418/474 


Abbreviations: CC, case-control study; TDT, tra 


ismlsslon/disequ 


librlum test. 









These findings are biologically plausible on some level. rs2254546 
is located at an intergenic region between FAM167A and BLK at 
8p22-23, a region harboring SNPs associated with several auto- 
immune diseases, such as systemic lupus erythematosus^*"'^', rheum- 
atoid arthritis^" and systemic sclerosis^"'™. PAM167A encodes a 
ubiquitously expressed gene but has not been functionally character- 
ized'". BLK encodes B-lymphoid tyrosine kinase, a Scr family tyro- 
sine kinase expressed primarily in the B-cell lineage, and transduces 
signal downstream following stimulation of B-cell receptors^*""". 
Little is known about the function of rs2254546 so far, and all that 
we know are the statistical associations of this polymorphism with 
j^p2o,23 p]-ji^j]-y Sjogren's syndrome'"'. However, Onouchi et al.^" 
have demonstrated that two other SNPs, rsl3277113 and rs2736340. 



which have been repeatedly associated with other autoimmune dis- 
eases, were in high linkage disequilibrium with rs2254546 in 
Japanese (D' = 0.99, r^ = 0.87; D' = 0.99, r' = 0.91). It has been 
indicated that the risk A allele of rs 13277 11 3 was associated with 
reduced expression of BLK mRNA in transformed B-cell lines^*". 
Meanwhile, Chang et al. have also provided evidence that the risk 
T allele of rs2736340 was associated with lower expression of BLK in 
peripheral blood B cells during the acute stage of KD, and decreased 
BLK expression in peripheral blood B cells may alter B cell function 
and predispose individuals to KD". AU of the above indicate the 
reasonability of our results. What's more, it is well known that most 
of the variants identified by GWASs cannot be per se causal but imply 
the probability of being in linkage with the "real" causal variants''', so 



Inverse Variance (Fixed-Effects) ORs 



Dnouctii (GWAS), 2012 
Onouclii (CC1), 2012 
Onouchi {CC2), 2012 
Onouchi (TDT1), 2012 
Onouchi (TDT2), 2012 
Onouchi (TDT3), 2012 
Onouchi (TDT4), 2012 
Yan, 2013 
Current study, 2013 
Pooled OR 




1 \ — \ — r 

2.0 2.5 3.0 4.0 



OR( 95 % CI) 

Figure 1 | Forest plot of the association between rs2254546 G allele and KD risk. For each study, the estimate of OR and its 95% CI is plotted with a box 
and a horizontal line. Fixed-effects pooled OR = 1.55, 95% CI = 1.42-1.70, P < 0.001; Q = 13.00, Pheterogeneity = 0.112. 
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Table 3 | Meta -analysis of rs2254546 in association with KD risk stratified by study design 



Case-control TDT OR (95%CI) 



Study ID 


Cases G/A 


Controls G/A 


Transmitted G 


Non-transmitted G 


Case-control 


TDT 


Onouchi (GWAS), 2012 


702/154 


4866/1892 


- 


- 


1.77(1.48-2.13) 


- 


Onouchi (CCl), 2012 


761/179 


522/234 






1.91 (1.52-2.38) 




Onouchi (CC2), 2012 


454/1 14 


808/330 






1.63 (1.28-2.07) 




Onouchi (TDTl), 2012 






75 


60 




1 .25 (0.89-1 .76) 


Onouchi (TDT2), 2012 






24 


16 




1.50 (0.80-2.82) 


Onouchi (TDT3), 2012 






21 


12 




1.75 (0.86-3.56) 


Onouchi (TDT4), 2012 






25 


19 




1.32 (0.72-2.39) 


Yan, 2013 


590/126 


1246/384 






1.44(1.15-1.81) 




Current study, 201 3 


664/172 


724/224 






1.19(0.95-1.50) 




Pooled 


3171/745 


8 1 66/3064 


145 


107 


1.57(1. 34-1 .85)t 


1.35 (1.05-1.74)* 



'Random-effects pooled OR, P < 0.001 ; O = 1 0.84, Pheterogeneiiy = 0.028 
'Fixed-effects pooled OR, P = 0.01 8; Q = 0.83, Ph,„,„,og,„,iiy = 0.843. 



we can speculate that rs2254546, like the two relevant SNPs, 
rs2736340 and rsl3277113, is possibly the real causal variant, but 
more likely to be just in linkage equilibrium with the real causal 
variant. 

The current study had its own special advantages. Firstly, we 
adopted a comprehensive analysis strategy (a case-control study 
and followed by a meta-analysis) to discern the association between 
rs2254546 and KD susceptibility; Secondly, in this meta-analysis, 
significant associations were detected in both case-control and 
TDT subgroups, which confirmed the statistical consistency of the 
two study designs. The study power was strongly enhanced by integ- 
rating these two study designs together due to the enlargement of 
sample size. Thirdly, the between-study heterogeneity in the overall 
meta-analysis was so negligible that these studies were considered 
compatible. Finally, the results of cumulative analysis, sensitivity 
analysis and publication bias further provided statistical evidence 
of the high quality of the meta-analysis. 

Despite of the clear advantages we have mentioned above, several 
limitations should be figured out at the same time. Firstly, our own 
case-control study failed to replicate the significant association in 
allelic model, and the relatively small sample size and low statistical 
power might be taken into account. Secondly, the association 
between rs2254546 and KD risk is merely on statistical level, so the 
lack of functional experiments makes its role in KD pathogenesis still 
an enigma. Thirdly, KD is a complex disease caused by both genetic 
and environmental factors. Fiowever, the insufficient information 
about the environment limited the further evaluation of gene-envir- 
onment interaction. 

In summary, our case-control study combined with the following 
meta-analysis might provide convincing evidence for the significant 
association between rs2254546 polymorphism and KD susceptibility. 
Fine-mapping and sequencing studies, as well as functional studies are 
further warranted to identify and validate the potential causal variants. 



Methods 

Study participants. Tills study contained 428 children diagnosed witli KD, all of 
whicli were unrelated ethnic Han Chinese and were consecutively recruited between 
April 2009 and September 2012 from Children's Hospital, Zhejiang University School 
of Medicine, China. The diagnosis of KD was based on the 5th revised edition of the 
guidelines established by the Kawasaki Disease Research Committee in Japan in 2002 
(http://kawasakl-dlsease.org/diagnostlc/). The controls consisted of 478 gender- 
matched unrelated healthy Han Chinese children without any evidence of infection, 
part of which have ever been included in our previous study"*^. AU the controls were 
recruited from children who took the entrance physical examination for primary 
school in Children's Hospital, Zhejiang University School of Medicine. 

This study was approved by the ethics committee of the Children's Hospital, 
Zhejiang University School of Medicine, and the methods were carried out in 
accordance with the approved guidelines. Informed consent was obtained from all 
participants or parents/careglvers of aU the subjects who were studied. 

Genotyping. Genomic DNA was extracted from 2 ml peripheral blood sample that 
was collected from each participant, using the RelaxGene Blood System DP319-02 
(Tlangen, Beijing, China) according to the manufacturer's protocol. The 
concentration and the optical density of DNA was confirmed by NanoDrop 1000 
spectrophotometer {Thermo Fisher Scientific, Waltham, Massachusetts, USA). 
Genotypes were determined by TaqMan SNP Genotyping Assay (Applied 
Biosystems, Foster City, CA, USA) as described below, which was performed without 
knowledge of participants' disease status. The final reaction volume per well was 5 \iL 
for a 384-weU plate, including 20 X working stock of SNP Genotyping Assay 
(0.25 nL), TaqMan Universal PGR Master Mix (2X), No AmpErase UNG (2.5 piL), 
DNase-free water (1.25 jiL), and Template DNA (1 [iL). The first two of which were 
customized by Applied Biosystems. The reaction was carried out with an initial step of 
10 minutes at 95 'C, which was the optimal temperature for AmpllTaq Gold Enzyme 
Activation; then followed by 50 cycles of 15 seconds at 92' C and 90 seconds at 60 'C. 
After PGR amplification, allelic discrimination plate read and analysis were 
performed using Applied Biosystems 7900HT Fast Real-Time PGR System, followed 
the manufacturer's protocol. 

For quality control, a 5% random samples were tested twice to assess the repro- 
ducibility, with a concordance rate of 100%. A total of 892 samples were genotyped 
successfully, so that the genotyping call rate of rs2254546 was 98.45%. 

Statistical analysis. HWE was assessed by the goodness-of -fit test for genotypes in 
the control group. Pearson's test was applied to estimate the differences in gender 
and distribution of genotypes between cases and controls. After adjusting for gender, 
unconditional multivariate logistic regression was employed to estimate ORs and 



Table 4 | Meta-analysis of rs2254546 in association with KD risk in Asian population 



Study ID 


Case-control 




TDT 


OR (95%CI) 


Cases G/A 


Controls G/A 


Transmitted G 


Non-transmitted G 


Onouchi (GWAS), 2012 


702/154 


4866/1892 






1.77(1.48-2.13) 


Onouchi (CCl), 2012 


761/179 


522/234 






1.91 (1.52-2.38) 


Onouchi (CC2), 2012 


454/1 14 


808/330 






1.63 (1.28-2.07) 


Onouchi (TDT2), 2012 






24 


16 


1 .50 (0.80-2.82) 


Yon, 2013 


590/126 


1246/384 






1.44(1.15-1.81) 


Current study, 2013 


664/172 


724/224 






1.19 (0.95-1.50) 


Pooled 


3171/745 


8 1 66/3064 


24 


16 


1.58 (1.44-1.74)t 


'Random-effects pooled OR, P < 0.001 ; O = 


10.87, n,e,e,oae„eiV = 0.054. 
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Table 5 j Sensitivity analysis of allelic model 


Study omitted 


Q 


^heterogeneity 


OR (95%CI) 


P 


Dnnnrh'i tC^\AJA's\ 9019 


1 n ?9 

1 VJ.OZ 


0 171 
VJ. 1 / 1 


1.49 ( 


1.34-1.65) 


•^Ci 00 1 

^W.KJKJ 1 


Onouchi (CCl), 2012 


9.17 


0.241 


1.49 ( 


1.35-1.65) 


<0.001 


Onouchi (CC2), 2012 


12.83 


0.076 


1.51 ( 


1.31-1.75) 


<0.001 


Onouchi (TDTl), 2012 


1 1.31 


0.126 


1.58 ( 


1.44-1.73) 


<0.001 


Onouchi (TDT2), 2012 


12.98 


0.072 


1.53 ( 


1.34-1.75) 


<0.001 


Onouchi (TDT3), 2012 


12.89 


0.075 


1.52 ( 


1.34-1.74) 


<0.001 


Onouchi (TDT4), 2012 


12.69 


0.080 


1.54 ( 


1.35-1.76) 


<0.001 


Yan, 2013 


12.51 


0.085 


1.54 ( 


1.33-1.79) 


<0.001 


Current study, 20 1 3 


6.76 


0.454 


1.63 ( 


1.48-1.80) 


<0.001 



their 95% CIs for the effect of the SNP on KD susceptibility. To avoid the assumption 
of genetic models, dominant, recessive and additive model were all analyzed 
respectively. All above statistical analyses were performed in SPSS V13.0 and all P 
values were two-tailed with a significant level of 0.05. 

Meta-analysis of rs2254546 in association with KD risk. To further confirm the 
association between rs2254546 and KD susceptibility, a meta-analysis combining 
published literatures and our current study was conducted. We searched various 
databases including PubMed, EMBASE, ISI Web of Science and Google Scholar up to 
October 2013 without language restriction, using the keywords 'rs2254546 or 
FAM167A-BLIC combined with 'Kawasaki disease'. We further checked five large 
Chinese biomedical bibliographic databases to expand the coverage of our searches^^. 
References listed in the retrieved articles and reviews were also scanned for missing 
information. The included articles had to meet the following criteria: (1) either case- 
control study or TDT design; (2) availability of allelic frequency necessary for 
calculating ORs and their 95% CIs; (3) clear definition of KD, either based on the 
diagnostic criteria from the Committee on Rheumatic Fever, Endocarditis, and 



Kawasaki Disease, Council on Cardiovascular Disease in the Young, American Heart 
Association^, or the diagnostic guidelines for Kawasaki disease (the 5^' revised edition, 
http://kawasaki-disease.org/diagnostic/}, published by the Kawasaki Disease 
Research Committee of Japan; (4) studies on humans. When multiple articles had 
overlapped study participants, the one with more complete design or larger sample 
size was fmally selected. 

For each of the eligible studies, the following information was extracted indepen- 
dently by two reviewers {W.W. and J.L.): first author's name, year of publication, 
country, ethnicity of study participants, study design, study type, genotyping method, 
sample size, frequencies of alleles in case and control group in case-control studies 
and numbers of transmitted and non-transmitted alleles from heterozygous parents 
to affected offspring in family-based studies. The Cochran's /^-based Q statistic test 
was applied to assess between-study heterogeneity. As regard of the p-value cut-off, 
both the use of a significance level of 0. 1 or 0.05 are commonly accepted in this type of 
|.gg|.37,38 jj^jg study, we chose 0.1 as the significance p-value cut-off for the het- 
erogeneity. If P < 0.1, a random -effects model, using DerSimonian and Laird 
method^^, was employed to compute the pooled ORs; otherwise, a fixed-effects model, 



Cumulative Meta-Analysis: 
Inverse Variance (Fixed-Effects) ORs 



study Added 
Dnouchi (GWAS), 2012 

Study Added 
Onouchi (001), 2012 

Study Added 
Onouchi (002), 2012 

Study Added 
Onouchi (TDTl), 2012 

Study Added 
Onouchi (TDT2), 2012 

Study Added 
Onouchi (TDT3), 2012 

Study Added 
Onouchi (TDT4), 2012 

Study Added 
Yan, 2013 

Study Added 
Ourrent study, 2013 




0R( 95% CI) 

Figure 2 | Forest plot of cumulative meta-analysis between rs2254546 G allele and KD risk in chronological order. The estimate of ORs and their 95% 
CIs are plotted with boxes and horizontal lines. 
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-1.00 



-0.50 



0.00 



0.50 



LogOR 



1.00 



1.50 



Figure 3 | Funnel plot of natural logarithm of OR against inverse 
standard error in each study. Regression test for funnel plot asymmetry P 
= 0.815. 



using Mantel -Haenszel method'*'', was adopted. For the synthesis of case-control 
studies and TDTs, a method described by Kazeem and FarralP* was used to obtain a 
combined estimate of the magnitude of the association. The "Catmap" software*^ for 
this method was downloaded from the comprehensive R archive network http:// 
www.r-project.org. On condition that the heterogeneity existed, a meta-regression 
analysis would be conducted to explore the possible sources of the heterogeneity 
across studies^. Stratified analysis was carried out subsequently, on the basis of the 
suspected sources of heterogeneity reported by the meta-regression analysis. 
Afterwards, sensitivity analysis was performed to assess the influence of each study on 
the overall estimate. Cumulative meta-analysis was also conducted by adding studies 
one at a time in chronological order. Funnel plot and regression test were utilized to 
assess publication bias. All statistical analyses and figures were performed in R soft- 
ware V2.13.2 and P values less than 0.05 were considered statistically significant for all 
tests except for Q test for heterogeneity. 
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